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A certain type of integer grid, called here an echelon grid, is an 
-sj , object found both in coherent systems whose components have a fi- 

nite or countable number of levels and in algebraic geometry. If a = 
(ai,...,Qd) is an integer vector representing the state of a system, 
pH , then the corresponding algebraic object is a monomial x"^ ■ ■ ■ x'^'^ in 

^0 ' the indeterminates x\,. . . , Xd- The idea is to relate a coherent system 

to monomial ideals, so that the so-called Scarf complex of the mono- 
mial ideal yields an inclusion-exclusion identity for the probability 
of failure, which uses many fewer terms than the classical identity. 
Moreover in the "general position" case we obtain via the Scarf com- 
plex the tube bounds given by Naiman and Wynn [J. Inequal. Pure 
Appl. Math. (2001) 2 1-16]. Examples are given for the binary case 
's^ . but the full utility is for general multistate coherent systems and a 

>^^^ ' comprehensive example is given. 

vQ ■ 1. Introduction. The study of network reliability has received increasing 

(^ I attention in the recent decades because of its applications to computer net- 

'nT ■ works and communication systems. Shier [21] points out that one of the most 

■5^ . interesting aspects of this subject is the "variety of discrete, combinatorial, 

1 -^ I and algebraic mathematics that can be found lurking just underneath this 

C^ ' practical veneer. . . . The new approaches developed in the context of network 

reliability can have ramifications beyond that particular venue." 
^ I Naiman and Wynn [16] established a connection between reliability struc- 

tures and some special topological objects, called abstract tubes. This study 
led to a second paper [10], where the interactions with computational alge- 
C^ ' bra and Grobner bases were investigated. This was motivated by the work 

of Diaconis and Sturmfels [4] that initiated the use of algebraic techniques 
in probability and statistics. A section on reliability was also included in the 
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monograph ([17], Section 4.4) devoted to establishing or deepening the hnks 
existing between algebra and statistics. The work of Dohnien (see, among 
others, [5, 6, 7]) concentrates more on the combinatorial approach but fruit- 
ful and interesting synergies with the abstract tube theory and the algebraic 
approach are foreseen. 

In the present paper the authors investigate a strong link between reliabil- 
ity theory and certain algebraic structures which we believe will eventually 
unify much of the work just mentioned. The two particular objects covered 
in this paper are the Scarf complex and the finely graded Hilbert series. 
These can contribute to finding more efficient computations and identities, 
similar to those in [16] and [10]. 

The power of the suggested approach lies in the fact it can be gener- 
ally applied to any kind of coherent multistate system without any further 
restrictions. Former approaches proposed in the literature are limited to spe- 
cial cases such as binary systems or source-to-terminal reliability for planar 
networks. Moreover, because of the link with the tube theory the method 
immediately gives bounds which are only available, outside the standard 
Bonferroni case, for certain such cases. Section 5.2 contains a brief discus- 
sion of the existing literature and provides examples to illustrate the new 
methodology. 

2. An overview. The main point of contact between the disciplines of 
algebra and reliability theory is an echelon grid of points with nonnegative 
and integer coordinates. Such a grid (or lattice) can identify a coherent 
reliability system but also plays a specific role as an algebraic object. An 
example is given in Figure 1. 

A reliability system of d components is a system whose failure or non- 
failure state is determined by the state of each of its components. We deal 
here with the case where each component can assume a discrete and pos- 
sibly countable number of states {0, 1, 2,3,.. .} corresponding to increasing 
levels of efficiency. We refer to such systems as multistate systems. Section 
5.5 deals with the case of a continuous distribution on the component per- 
formance and more formal definitions about reliability systems are given 
and discussed in Section 5. A state of the given system is described by a 
nonnegative integer vector of length d. The state space D (P = N in this 
general setting) can therefore be decomposed into the set of configurations 
that corresponds to failure of the system and the set of configurations that 
corresponds to nonfailure. We indicate these two sets with JF and J^, respec- 
tively. Obviously T> = J-U J- (i.e., J- is the complement of JF in P). 

Coherent systems are systems for which improving the state of a single 
component cannot lead the system from a nonfailure (or operating) state 
to a failure state. Equivalently, degrading a component state cannot bring 
the system from failure to nonfailure. Thus, in the integer representation 
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a coherent system is such that if a state point a belongs to the failure set 
J-, then any other point with component state levels worse (less) than a 
must belong to J-. Similarly the nonfailure set, being the complement of J-, 
must have a similar property: if a is a nonfailure point, then any point with 
coordinates greater than or equal to a's coordinates must still be a nonfailure 
point. This is equivalent to saying that the integer representation of the 
failure set must have an echelon structure; that is, the grid representing J^ 
(or equivalently !F) must have no "holes." 

For example Figure 1 represents a two-component coherent system where 
the hollow dots correspond to failure states and the filled dots correspond 
to operating (nonfailure) configurations. 

The link with algebra is the construction of the monomial x" = x"^X2^ ' ' ' ^d 
with the point a = {ai,a2, ■ ■ ■ , ctd) given by its exponent vector. For exam- 
ple, the filled dots in Figure 1 represent the monomials belonging to a certain 
monomial ideal M. For the introductory material in algebra we have mainly 
used [3, 12, 14] as references. The echelon property, equivalent to coherence, 
translates into the algebraic language as the "order ideal" property: if the 
monomial x° belongs to a monomial ideal, then so do all the monomials x^ 
such that x" divides x^, that is, such that a ^ /? : aj < Pi,i = I,. . . ,d. The 
corner points of the set T (those labeled with their coordinates in Figure 
1) are the minimal generators of this monomial ideal but it is interesting to 
point out that these points play also a specific role as reliability theory ob- 
jects. They are called minimal nonfailure points in Section 5. In this paper 
we shall only study coherent systems. 



"d 
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Fig. 1. Monomials in M = {x\,x\x\,x'2) (black dots) and monomials not in M (hollow 
dots). 
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The main observation, which we will develop in the following sections, is 
that there are some points on the boundary between the failure set and the 
nonfailure set that are fundamental to the identification of the points in T 
and in T . We will see that these special points form a simplicial complex, 
called the Scarf complex, which then plays a key role in both the algebraic 
geometry and reliability theory. 

Sections 3 and 4 give an introduction to the algebraic setting and Section 

5 establishes and formalizes the connection with reliability. 

3. Monomial ideals. 

3.1. Some basics. Let IC be a field and let S = M.[x] = ]K[xi, . . . , Xd] be the 
polynomial ring in d indeterminates. A monomial in M.[x] is a product x" = 
x"^X2^ • • • x^'' for a vector a = (ai, . . . , a^) € N'' of nonnegative integers. 

A monomial ideal M is an ideal of S generated by some monomials x", 
for a belonging to a subset ^ C N'^ (possibly infinite). The ideal M consists 
therefore of all the polynomials of the form J2aeA haX°^, with ha € S. In this 
case we write M = (x"|a G A). 

The following lemma allows us to characterize all the monomials that lie 
in a given monomial ideal. 

Lemma 1. Let M = (x"|a ^ A) he a m,onomial ideal. Then a m,onom,ial 
x^ lies in M if and only if x^ is divisible by x" for some aG A. 

Observe that the set 

a + N"* = {a + 7 I 7 G N"*} 

consists of the exponents of all monomials divisible by x". This observation 
allows us to draw pictures of the monomials in a given monomial ideal and 
visualize this set as a union of positive and integer coordinate points. For 
example, the filled dots in Figure 1 represent all the monomials in the ideal 
M={xl,xlxl,xl). 

A polynomial / is in a monomial ideal M = {x" \ a (z A) if and only if 
each term of / is divisible by one of the given generators x". From this it 
follows that a monomial ideal is uniquely determined by its monomials; that 
is, two monomial ideals are the same if and only if they contain the same 
monomials. For details and proofs see [3, 14]. 

The main result in this section is given by the following theorem (Dick- 
son's lemma), which states that each monomial ideal is uniquely and finitely 
generated. 

Theorem 1. A monomial ideal M = (x°^|a S ^) C K[xi, . . . ,Xd] can be 
written in the form M = (x"^, . . . ,x°^') , where Oi G A. In particular, M has 
a finite basis. 



SCARF COMPLEX FOR COHERENT SYSTEMS 5 

The polynomial ring S = ]K[x] can be seen as a K-vector space and there- 
fore it can be decomposed into the direct sum S = ^aeN^- '^a, where Sa is the 
K-span of the monomial x*^. Since Sa ■ Sf^ Q Sa+p, we say that S is an N*^- 
graded IK-algebra. More generally, an S-module M is said to be Z'^-graded 
if M = 0^g2d Mp is a direct sum of K- vector spaces with Sa • Mp C Ma+p ■ 
For definitions and properties of modules of polynomial rings see [8] . 

A monomial ideal M and the corresponding factor ring S/M are both 
Z'^-graded 5- modules: M = ^^a^j^ Sa and S/M = ^^a^j^.f Sa- 

The hollow dots in Figure 1 within the nonshaded area form a K-basis 
for S/M . We show in the next section that the Hilbert series H{S/M; x) is 
the sum of all the monomials that are not contained in M. 

Given a Z^'-graded module M, the Z^'-graded shift M[a\ for a € Z'^ is 
the Z'^-graded module defined by M[a]^ = Ma+p- In particular, the free S- 
module of rank 1 generated in degree a is S'[— a]. There is an isomorphism 
between S'f— a] and the principal ideal {x") C S. 

3.2. Hilbert series and free resolutions. Given a Z'^-graded module Ma 
we consider its dimension as a vector space on IK and we call it dini^{Ma)- If 
diin^{Ma) is finite for all a G Z*^ the (finely graded) Hilbert series is defined 
as the formal power series 

H{M;x):=HiM;xu...,Xd)= J2 dimK(M,)x". 



For example, we have that 

HiS;x)= E^" = nY 



d ^ 



(namely the sum of all monomials in S), and for a G Z'^, 
H{S[-a]]x)= Y. ^^ 

= x°(l + xi + xl + ■■■)■■■ {I + Xd + xl + ■■ ■) 

1 — Xl 1 — Xrf 

x" 



nf=i(i-^.) 

For a monomial ideal M we have that 

H{S/M;x)= Y. ^"' 

qGN<*\M 

namely the sum of all monomials not in M. 
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A homological complex of S-modules is a sequence 

of S'-module homomorphisms such that (pi-i o (^j = 0. A complex is exact at 
the ith. step if it has no homology there, that is, if Kernel((^j_i) = Image(i;^j). 
The complex is exact if it is exact at the ith step for all i £Z. 
A free resolution of an S'-module M is a complex 

of free modules which is exact everywhere except the 0th step, and such that 
M = Coker((/)i) = Fo/Image((/)i). 

Every S-module has a free resolution, with length less than or equal to 
d. If M is Z'^-graded, then it has a Z'^-graded free resolution. If, in addition, 
M is finitely generated, there is a Z'^-graded resolution M <— F in which 
all the ranks of the Fi are finite and simultaneously minimized. Such an F 
is called a minimal free resolution of M and is unique up to noncanonical 
homomorphism. 

Given a short exact sequence <— M" <— M <— M' ^- 0, the rank nullity 
theorem of linear algebra implies that dimK(Ma) = dimK(M^') + dimK(M^) 
for all a, and hence H{M;x) = H{M";x) + H{M';x). More generally, if 

^ M ^ Fo ^ Fi < 

is a finite sequence such as a free resolution, then 

(2) H{M;x) = Y,{-iyH{Ff,x). 

i 

In particular, if M is finitely generated, the existence of a finite-rank free 
resolution for M implies that the Hilbert series of M is a rational function 
of X, because it is an alternating sum of Hilbert series of S'[— a] for various 
a. Moreover the denominator can always be taken to be Yli{l — Xi). 

We show, with an example, the connection between minimal resolution 
and Hilbert series and why we are particularly interested in it. 

Example 1 . Consider a monomial ideal in two variables: M = {x""^ y^ , x""'- 
y*2 , . . . , x°"'y^'~) in S = ]K[x, y] . As mentioned above the finely graded Hilbert 
series of the factor ring S/M gives the sum of all the monomials not in M 
(and for the monomial ideal M it gives the sum of all the monomials in M) . A 
way to obtain the Hilbert series for S/Ad is to proceed by inclusion-exclusion, 
subtracting from the sum of all the monomials in S the list of all the mono- 
mials in the quadrants (x'^'y**) for all i G {1, . . . ,r}, then adding back the 
monomials in the two-way "intersections" {lcm.{x°'*y^%x°'^y'^^)) = {x°'*y^^) H 
(x"^y^^), then removing the three-way intersections (x°'y*') n {x°'^y^^) n 
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(x"''' y^* ) , and so on. In this way, after r steps we have counted ah the mono- 
mials the right number of times. However, for a large number of generators 
this procedure is far from efficient since many terms cancel out. It is easy 
to see that in this two-dimensional case it is enough to add back only the 
principal ideals obtained from couples of adjacent generators. 

In the literature the highly nonminimal free resolution of S/M given by 
the inclusion-exclusion process is called the Taylor resolution, which we de- 
scribe in more general terms in Section 4.1. The most compact representation 
yields instead the minimal resolution. Unfortunately a construction of the 
minimal free resolution has not been found for all arbitrary monomial ideals. 
In the literature there are only two general constructions for resolving arbi- 
trary monomial ideals: the Taylor resolution and Lyubeznik's subcomplex. 
We show in Section 4, based on [2, 13, 14], a construction of a minimal free 
resolution for a specific class of monomial ideals, called generic monomial 
ideals. Also a nonminimal free resolution of S/M based on deformation of 
the exponents is obtained for nongeneric monomial ideals. Even though this 
resolution is not minimal it is generally much smaller than Taylor's. We 
show in Section 5 that in reliability theory the generic monomial ideal has 
a very natural interpretation but even for binary systems (which typically 
lead to the nongeneric case) the deformation procedure still gives excellent 
results. 

Notice that the deformation procedure is gratifyingly similar to the per- 
turbation presented in [16] and [10], and the nonminimal resolution obtained 
via the Scarf complex after deformation leads to the abstract tube formula 
described in [16] and [10]. 

4. Monomial resolutions. In this section we first describe the construc- 
tion of a (highly nonminimal) resolution for an arbitrary monomial ideal, 
called Taylor resolution. Then a method is described to obtain the minimal 
resolution for generic monomial ideals, and finally the deformation procedure 
to deal with arbitrary monomial ideals. For more details, see [2, 13, 14]. 

4.1. Taylor resolution. For a monomial ideal M = {mi, . . . ,mr) and a 
subset / C {1, . . . ,r} we set m/ = lcm(?7ii|i € /). Let q/ G N be the expo- 
nent vector of mj and let S{—ai) be the free S-module with one generator 
in multidegree aj . Let Fg be the K- vector space whose basis elements e/ cor- 
respond to the index sets / C {1, . . . ,r} of length s. Define the differential 
ds-.Fs^ Fs^i by 

ds{ei) = J2^^ST^i^^^) e/V' 

where sign(z,/) is (—1)-''^^ if i is the jth element in the ordering of /. It is 
possible to show (see [8]) that 

O^Fo^Fii ^ Fr-l^Fr ^ 
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is a free resolution of S/M of length r, and it is called Taylor resolu- 
tion of S/M . The resolution can be also seen as the Z'^-graded module 
F = 0/c{i r\^{~^i) that has 2^ terms and it is therefore very far from 
minimal if r ^ n. 

It is natural to write the Hilbert series for a Taylor resolution of the factor 
ring S/M as follows: 

H{S/M;x,y) 

= H{S;x,y)-H{M;x,y) 

= E^"- E ^"+ E x° + ... + (-ir E ^^ 

ogN<* x°' e (mj) x" G (m/> x°'e{nij) 

IC{l,...,r} IC{l,...,r} I={l,...,r} 

\I\ = l |/| = 2 

(3) 

where (m/) is the principal ideal generated by the monomial mi. Similarly 
for the monomial ideal M, 

H{M-x,y) 

= E ^"- E x° + ... + (-1)^+1 Y. ^"- 

^ ' x°' e (mi) x" e [mi) x°' G (m/) 

/C{l,...,r} /C{l,...,r} 7 = {l,...,r} 

|/| = 1 |/|=2 

Notice the analogy between the structure of the above formula and the 
classical inclusion-exclusion formula used in reliability. Simply replace any 
term x" by the indicator function of the orthants Qaj = {/? | /? > a/}- 

We show in the next sections that the Scarf complex allows a more efficient 
formulation of the above expression and yields exactly the tube formula 
proposed in [16] and in [10]. For clarity we return to the example of Figure 
1. 

Example 2. We construct now the Taylor resolution and the Hilbert 
series of S/M where M = {x^,x'^y'^, y^) is the monomial ideal given in Section 
3 and represented in Figure 1. 

The required resolution is given by 

O^Fo^Fi^F2^F3^0. 

The dimension of each vector space Fg {s = 0, 1, 2, 3) in the sequence is given 
by (^) (i.e., the number of subsets I C {1, . . . , d} of length s). 

The differentials dg can be expressed via monomial matrices where each 
column corresponds to an index set / and contains the vector ds{ei). There- 
fore we obtain the following: 



Q^. [ ^3 ^2„2 „3 



X x y y^ 

/={!} /={2} /={3} 
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d.: 



x^y"^ 







y 







y 



7={1,2} -r={l,3} -r={2,3} 



2-3^2 

/={1,2,3} 

3 3 

mj^x y 

For example, for / = {1, 2} it is 

mi2 



92(612) 



mi 



-ei 



^12 

7712 



-62 







which corresponds to the first column of the matrix representing 82- 
Now from (1) and (2) it follows that 



H{S/M;x,y) 



1 



x'^y'^ — y^ + x^y"^ + x^y^ + x'^y^ 



x y 



(l-x)(l-y) 



which is exactly the sum of monomials not in M, as expected. 

Notice that before cancellation of the term x^y^ the numerator of the 
Hilbert series for the Taylor resolution contains 2^ = 8 terms. 

4.2. Monomial resolution over labeled siniplicial complexes and the Scarf 
complex. It is possible to define a monomial resolution associated with a 
simplicial complex labeled by the generators of a given monomial ideal. 
In this section we show that the resolution of a generic monomial ideal 
associated with a specific simplicial complex, called the Scarf complex, is 
minimally free. The Scarf complex owes its name to Herbert Scarf, who 
introduced a similar complex in mathematical economics (see [20]). The 
importance of this structure for the resolution of a monomial ideal was only 
recently understood and explained in [2]. 

Definition 1. Let V = {vi, . . . ,Vr} be a finite set. A (finite) simplicial 
complex A on y is a collection of subsets of V such that Ii € A whenever /i C 
I2 for some I2 G A, and such that {vi} G A for i = 1, . . . , r. The elements of A 
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are called faces, and the dimension dim(/) of a face / is the number |/| — 1. 
The dimension of the simplicial complex A is dim(A) = niax{dini(/) : / G 
A}. 

Note that the empty set is a face of dimension —1 of any nonempty 
simplicial complex. Faces of dimension and 1 are called vertices and edges, 
respectively. The maximal faces under inclusion are called facets of the sim- 
plicial complex. 

For a monomial ideal M in 5" = ]K[xi, . . . ,Xd] we can construct a special 
simplicial complex whose vertices are labeled with the generators of M. 
Additionally we label each face / of the simplicial complex by the least 
common multiple nij = lcm(77ij|i € /) of its vertices and we restrict to the 
faces with unique labels. This leads to the following definition. 

Definition 2. For a monomial ideal M = {nii, . . . , m,r) we call the Scarf 
complex the simplicial complex Am consisting of sets of minimal generators 
with unique labels 

Am = {I Q {1, . . . ,r} \ nij ^ nij for all J Q {1, ■ ■ ■ ,r} other than /}. 

It is not difficult to see that a Scarf complex in S" = IC[xi, . . . , Xd] is still a 
simplicial complex and has dimension at most d— 1. 

For i € Z let Fi{AM) be the set of i-dimensional faces of Am and let 
^Fi(AM) |-jg g^ IK-vector space whose basis elements ej correspond to the i- 
faces / G Fj(AAf). Let Fa^^ be the so-called N'^-graded chain complex of Am 
over S obtained as 

^ IK^-i(^m) ^...^ ]^F,^i{^m) d^-^F,{AM) , ^jzi ]K^n-l(AM) ^ Q 

with differentials as in Section 4.1, 

5s(e/) = Vsign(i,/) ei\i. 

For a special class of monomial ideals called generic, the complex defined 
by the simplicial complex Am is minimally free as Theorem 2 states (for the 
proof see [2]). 

Definition 3. A monomial ideal M is called generic if no variable Xi 
appears with the same nonzero exponents in two distinct minimal generators 
of M. 

Note that in [13, 14] the authors prove that it is possible to define gener- 
icity in a weaker way and still achieve the same results that we are going 
to present. They therefore refer to the definition we give above as "strong 
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genericity." For the purpose of this article the stronger version is not only 
sufficient but also easier and more appropriate to use for application to relia- 
bility. The requirement that the monomial ideal is generic might seem quite 
strong. In practice though almost all monomial ideals are generic, in the 
sense that those which fail to be generic lie on finitely many hyperplanes in 
the space of exponents. From the reliability point of view the requirement is 
not too strong in the continuous distribution case (see Section 5.5), while in 
the binary case the deformation procedure can be used to obtain a resolution 
(typically nonminimal) as we shall see from examples. 

Theorem 2. For a generic monomial ideal M the complex Fa^v^ defined 
by the Scarf complex Ajv/ is the minimal free resolution of S/M over S. The 
fi -graded Hilbert series of S/M {i.e., the sum of all monomials not in M) 
is 



(5) 



E/eAM(-l)''^' ^i 



(1 -xi)---(l -Xd) 
and there are no cancellations in the alternating sum in the numerator. 

Theorem 2 allows us to write the Hilbert series for the factor ring S/M 
(and for the monomial ideal M) in a much more parsimonious way than in 
the Taylor resolution. The best way to see the difference is to compare (3) 
and (4) with the following expressions for H(M;x,y) and H{S/M;x,y): 



H{S/M;x,y) = 


:E-"- E -" 




aeN'' x" e (mi) 




le^M 




\I\ = 1 




+ > ; x" + • • • + 1 




X- G {mi) 




le^M 




|/|=2 


H{M-x,y) = 


■ E -"- E 




x"' G (mj) x" e (mj) 




IgAm I&^m 




|/| = 1 |7| = 2 



(6) 

E -"; 

" e {mi) 
{l,...,r} 

(7) H{M-x,y)= ^ x"- ^ :," + ... + (_i)'-+i ^ 

/ = {!,. ..,r} 

The Taylor resolution can be obtained from the full simplicial complex 
(with cardinality 2'') given by the set of all subsets of {1, . . . , r}, while in the 
Scarf resolution the summation is done only on the subsets belonging to Am- 
Let us see now with an example how we can obtain the minimal resolution via 
the Scarf complex for the monomial ideal considered in Example 2 (Figure 
1). The example shows a "by hand" construction of the Scarf complex while 
for most of the other examples the authors made continual use of available 
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functions in CoCoA. The algorithms behind these functions are related to 
the construction of the Hilbert series (see [14]) and seem reasonably efficient. 
Further research is needed to establish theoretical or experimental results 
on formal computational complexity or actual run time. We do not develop 
the computational aspects in this paper. 



Example 3. The Scarf complex for the monomial ideal M = (x'^, x^y^,y^) 
can be obtained in the following way: 



/: 
/: 
/: 
/: 
/: 
/: 
/: 



.3. 



2 2 

X y ■ 

3. 



{1} ^ mi 
{2}^mj 
{3} -^ mi = y 



{1,2} ^mi =x^y'^ 



:xV 



x^y^. 



{1,3} ^mi 
{2,3} ^mj 
{1,2,3} ^mj 

Thus the Scarf complex Am is given by 

Am = {{1},{2},{3},{1,2},{2,3}}. 

The index subsets that do not appear in Am are those corresponding to 
the monomial x^y^, and notice that this is exactly the term that cancels in 
Example 2. 

The resolution is given by 







Jx^xV^ c<3. 



y' 





— X 




y 

-X-2 



S^ 



0, 



and the Hilbert series (which can be read off directly from the labels of the 
Scarf complex) equals 



H{S/M; x) 



I 



^3 _ ^2y2 _ y3 _|_ ^3y2 _|_ ^2y3 



{l-x){l-y) 



Thus the number of terms in the numerator of the Hilbert series is six. The 
gain in comparison with the Taylor resolution, which returned eight terms, 
is not large for this small example but it becomes huge in high-dimensional 
problems. 
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4.3. Deformation of exponents. As pointed out in Section 4.2, when the 
generators of the monomial ideal are in generic position the Scarf complex 
leads to the minimal resolution. For arbitrary monomial ideals the construc- 
tion described in this section can be used to produce a (typically nonmini- 
mal) resolution via deformation of the exponent vectors of the generators. 
There is no guarantee that the resolution so obtained is minimal, but the au- 
thors' experience is that there is still a very great gain in resolution over the 
original Taylor complex. The use of the deformation procedure combined 
with the use of the CoCoA macros mentioned in Section 4.2 seems very 
powerful and certainly superior to complete enumeration followed by can- 
cellation. The example presented in Section 5.4 shows a typical result where, 
even after deformation, a considerable reduction in the number of terms in 
the inclusion-exclusion formula is achieved (from 511 to 31). Another ad- 
vantage of the approach presented here is that, even in the cases where 
minimality is not guaranteed, the method, unlike other methods described 
in the literature, still provides improved reliability bounds, as described in 
Section 5.2. 

For an arbitrary monomial ideal M = (mi , . . . , nir), let {oi = {an , . . . , aid) \ 1 < 
i <r} he the exponent vectors of the minimal generators of M. Choose vec- 
tors Ej = (sji, . . . ,£id) G K for I <i <r such that, for all i and s ^ t, the 
numbers Oj^ + £is and au + eu are distinct and Oj^ + £is < au + ^it implies 
ciis < ciit ■ Each vector £i defines a monomial x^' = x^'^ • • • x^*'' with real ex- 
ponents. We define the generic monomial ideal M^ in a polynomial ring with 
real exponents as follows: 

Mg = (mix^^,m-22;^^) ■ • ■ ,m-rX^''). 

We call Mg a generic deformation of M. Let A^,/^ be the Scarf complex 
of Mg. We now label the vertex of Am^ corresponding to mtx^^ with the 
original monomial ttij. Let Tg be the complex of S-modules defined by this 
labeling of Am . Then the following result can be obtained (see the proof in 

[2]). 

Theorem 3. The complex F^ is a free resolution of S/M over S. 

There is a simple way to deform M. This is the deformation given in [2] 
and is equivalent to the perturbation proposed in [10]. The method consists 
in picking an integer v > r and deforming M using £ij = i/v. 

5. Algebra and reliability. 

5.1. Coherent systems. We consider here systems as defined in [17] and 
[1] but with some slight differences in notation. 
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Definition 4. A system is a set S of d components. We code their 
increasing efficacy levels with the integers {0, 1,2,.. .}. 

An outcome is a nonnegative integer vector of length d describing the 
state of each component. We call V the set of all possible outcomes. 

A failure outcome is an outcome which leads to failure of the system S. 

The failure set T is the set of all failure outcomes. The nonfailure set is 
the complement oi T m.T>:T = T)\T. 

In most cases it is natural to assume that replacing a component by a 
component at a higher efficacy level will not lead to a deterioration of the 
system. Systems for which such an assumption is valid are called coherent 
systems. We consider only this type of system. For example, if the outcome 
(1, 2, 1, 3, 0) is an operating (nonfailure) outcome then also (1, 2, 2, 4, 0) must 
be an operating configuration. In terms of the usual partial order relation < 
in U^, for a coherent system we write that if q S .iT, then f3 £ T for all /3 such 
that a ^ P (i.e., ai < Pi, i = 1,. . . ,d). In the same way, in a failure configu- 
ration, replacing an operating component with a failed one cannot improve 
the system. Therefore, if (0,1,0,0,0) is a failure outcome then (0,0,0,0,0) 
must also be. This observation leads to the following definition. 

Definition 5. The minimal failure points are maximal points in the 
partial order relationship defined by :<: 

T* = {aeT\$f3eT:a^f3}. 

The use of the word "minimal" in the definition is justified by the fact 
that in the binary network literature the minimal failure points (sometimes 
called minimal cuts) are the minimal set of components whose failure ensures 
the failure of the system. The collection of all the minimal failure points is 
called a minimal failure set and is indicated with the symbol J-* . 

Similarly the minimal points (according to :<) in the nonfailure set are 
called minimal nonfailure points (paths in network theory): 

As mentioned in the Introduction it is easy to see how the failure set for a 
coherent system can be represented by an integer grid. Furthermore this grid 
has an echelon structure since the system is coherent. A similar observation, 
but with reverse inequalities, can be made for the nonfailure set. 

The parallel between a monomial ideal and a nonfailure set, and between 
minimal nonfailure points and minimal generators, is now fully established. 

Proposition 1. Given a system S of d components with nonfailure set 
T and minimal nonfailure set T* : 
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1. The minimal nonfailure points in T* , seen as the exponent vectors of 
monomials m lC[xi, . . . ,Xrf], are the (exponents of the) minimal genera- 
tors of a monomial ideal M . 

2. The points in T represent the (exponents of the) monomials belonging to 
the monomial ideal M generated by the minimal failure points. 

3. The points in T represent the (exponents of the) monomials belonging to 
the factor ring S/M . 

Thus we allow ourselves to use interchangeably the expressions "nonfailure 
set ^" and "monomial ideal generated by ^*," or "minimal failure points" 
and "minimal generators," even though, strictly speaking, we are dealing 
with monomials in one case and their exponents in the other one. 

The main interest in reliability theory is to give a measure of the per- 
formance of a system S by evaluating the probability of the failure (or 
nonfailure) of the system. 

Definition 6. For a system S the reliability function 7^(«S') is the prob- 
ability of the nonfailure set ^; the unreliability U(S) is the probability of 
the failure set J-. Clearly, we only need to concentrate on one of the two 
performance measures, since U(S) = 1 — TZ(S). 

In the literature it is possible to find many techniques to obtain the relia- 
bility of a system, but the most common problem encountered is the compu- 
tational effort required to evaluate this probability for large size problems. 

One of the methods proposed in the literature to calculate the reliability 
function in the case of a finite state space is due to Moore and Shannon [15] 
and is based on state space enumeration. Thus the reliability function for a 
network G with nonfailure set J^ can be written as 

7^(5) = 5]/^MPrM, 

where T> is the full set of network states and I^(io) takes the value 1 when 
io belongs to J^, otherwise. The approach is impractical because the space 
V has cardinality equal to m'^, if d is the number of components and m the 
number of levels that each component can assume. 

The reliability function can be formulated in terms of the minimal non- 
failure points. For a minimal nonfailure point a E J^* we indicate with Qa 
the event "the system is in one of the states /3, with /3 y a." This can also 
be described as the event that the system state is a point in the orthant that 
has the point a as corner point. 
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Notice that the orthants Qa correspond exactly to the set of monomials 
belonging to the principal ideals (x"^). Thus the nonfailure set can be written 
as a union of the orthants Qa based on the minimal nonfailure points a G T* , 

T= U Qa. 

The reliability function can therefore be obtained as the probability of a 
union of orthants. Since this union is not disjoint the classical approach to 
this problem is to use the inclusion-exclusion formula, 



7e(S)=Prob( U Q, 



= 5]Prob(g„)- Y. Prob(Q„nQ„0 + --- 

+ (-l)''+^Prob(Q^nQ„/n---), 
where r = \J-*\ is the cardinality of the minimal nonfailure set. 

5.2. The method. To obtain the improved version of the inclusion-exclusion 
formula for a system S with minimal nonfailure set T* = {ai, . . . ,«,.} we 
only need the Scarf complex A^/ associated with the monomial ideal M 
generated by the monomials x'^^ , ■ ■ ■ , x"'"' . The minimal resolution and the 
Hilbert series provide the background for our calculations. The Scarf com- 
plex shows directly which least common multiples are needed to be included 
in the identity. We decompose the complex in terms of the dimensions of its 
faces and we write 

A = {Aq, Ai, . . . , Ad_i} withd<n, 

where A^ is the set of faces of A of dimension s + 1 (thus Aq is the set of 
vertices, Ai is the set of edges, etc.). 

From (7) we obtain the improved version of the inclusion-exclusion for- 
mula, 

7^(5)=Prob( U qA 
(8) = E Prob(Q„,)- Y. Prob(g,J + --- 

+ (-l)'^"i Y Prob(Q,J, 

m/eAd_i 

where mi = lcm(?TT-j|i G /) is as defined in Section 4.2 and a/ is the exponent 
vector of the monomial mi. 
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This formula is equivalent to the improved inclusion-exclusion formula 
obtained in [16] (tube identity) for orthant arrangements. The authors prove 
there that truncating the formula at even and odd level leads to upper and 
lower bounds for the reliability, 

r+l r 

Y^i-iy Y. Prob(Q„,)<7^(5)<5:(-l)^■ J2 Prob(Q„,), 

j=0 m/gAj j=0 m,j£Aj 

(9) 

<r <d —l,r even. 

In [16] it is proved that such bounds are always at least as tight as the 
bounds from truncating the usual inclusion-exclusion lemma. It is of some 
interest that these inequalities have not, to the authors' knowledge, been 
established in the algebraic literature. 

5.3. The binary case. The methods of this paper were introduced to be 
applied to any general multistate coherent system and in the next section 
we give an example. We discuss first, however, the standard binary case on 
which there is much literature. 

This literature divides into three broad classes: cases derived from net- 
works, special nonnetwork cases (such as k-out-oi-n) and arbitrary binary 
examples. For network and some special nonnetwork cases there are iden- 
tities and bounds competitive with those presented in this paper. The best 
published results contain identities with the maximum amount of cancella- 
tion in the classical inclusion-exclusion identities. An example is the work 
by Satyanarayana and Prabhakar (see [18, 19]). For bounds the recent re- 
sults of Dohmen [5, 6, 7] also using tube identities are based on a somewhat 
different construction and are specialized to certain classes of networks. 

As mentioned, the immediate problem with using the Scarf resolution for 
binary systems is that the genericity condition only holds in trivial cases 
and the deformation procedure has to be used. The algebraic theory tells 
us that this does not guarantee a minimal resolution and, as pointed out 
in [10], different deformations give rise to different Scarf complexes. Thus, 
some deformations can be better than others, in the sense that they lead to 
a smaller Scarf complex and therefore to a more concise reliability identity. 
The authors are currently investigating the possibility of identifying the 
"optimal" deformation, that is, the deformation that leads to the minimal 
resolution. As the general theoretical minimal resolution problem appears 
to be unsolved in algebra an early theoretical solution is not expected, but 
fast computational methods look promising. 

Notice that in the binary case if we assume independence among the 
component failures and we indicate by pi the probability of nonfailure of 
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component i, the probability of an orthant Qaj associated with the index 
set I in the Scarf complex is easily obtained as 

(10) Prob(Q„J = n Pi- 

aii=l 

To summarize, the strategy for the binary case which derives from this 
paper is as follows: (i) describe the failure set; (ii) perturb the minimal failure 
points to general position; (iii) derive the Scarf complex; (iv) consider the 
labels (monomials) associated with the Scarf complex; (v) obtain the orthant 
probabilities; (vi) derive the final identities or bounds. 

Example 4 (Binary network). We illustrate our method by considering 
an eight-component binary system defined by the following nine minimal 
nonfailure points: 

J"* = {(1,0, 0,0, 0,1, 0,0), (1,0, 0,1, 0,0, 1,0), (0,1,0,1,0,1,0,0), 

(1,0,0, 1,1,0, 0,1), (0,1,0,0, 0,0, 1,0), (0,0, 1,1, 1,1, 0,0), 

(0,1, 0,0, 1,0, 0,1), (0,0, 1,0, 1,0, 1,0), (0,0, 1,0,0, 0,0,1)}. 

This example is taken from [6, 7]. The corresponding monomial ideal in 
]K[xi, . . . ,X8] is generated by the monomials corresponding to the minimal 
nonfailure points: 

M = {xiXQ,XlX4,X7, X2XiXQ,XlXiX^X^,X2Xj, X^XiX^XQ,X2X^X^,X'iX^Xj , X^Xg) ■ 

After perturbation (with £i = i/W) and ranking, we obtain the following 
generic monomial ideal: 

/ fy^ ry^ ry O^ T* 'y. /y" 'y ^y 'y fy ry^ ry^ ry^ /y^ fy ry^ ry^ ry^ ry ' ry ly fy ry 

fiS2fi84'^8 47 '^ ^72 ^28^762 

fy fy 'y /y> /y* ry ry ry /y 'y ry ,-, ry ry „ ry ry 'y ry ry ry ry ry ry 'y 'y 

X-iX^XriX AXr:.XaX^XQ^ XlX2.^3*^4"*'5"*'6 7 8, X^X^X'iX AXr:^Xp.X'7Xg^ 

XiXoXAX^XaX^XQ^ X\X2XoX4,Xr.XQXY^ XoXo) . 

The corresponding Scarf complex A ~ contains 103 elements, which are the 
index subsets of the following six facets: 

{{12479}, {13689}, {12579}, {12589}, {13579}, {13589}} , 

where 1, 2, . . . , 9 are the labels of the minimal generators of M. 

The reliability function is obtained following (8) and (10) and it is based 
on the probabilities of the orthants identified by the Scarf complex. For ex- 
ample, the index set I = {127} corresponds to the least common multiple of 
the first, second and seventh minimal generators: mj = X1X2X4X5XQX7X8 = 
lcm.{xiXQ,xiX4XY,X2X5Xs). In the independence case, using (10), the proba- 
bility of the corresponding orthant is expressed in terms of the probability 
of nonfailure of the single components, pi,i = l,. . . ,8, 

ai = (1,1,0,1,1,1,1,1), Proh{Qaj) =PlP2P4P5P6P7P8, 
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where a/ is the exponent of the monomial mj. 

Notice that, interestingly, the Scarf complex here coincides exactly with 
the tube simplicial complex obtained by Dohmen [6, 7] in the case of a 
binary planar network whose reliability function is defined as source-to- 
terminal reliability. However, the Scarf complex method can be applied to 
much more generally defined coherent systems, as for example the binary 
nonnetwork case or the multistate case. 

Example 5 (Binary nonnetworks). The recent literature contains sev- 
eral efficient algorithms to calculate the reliability of special classes of co- 
herent systems. For example, the method presented in [19] gives an efficient 
algorithm that can be used to evaluate the source-to-terminal reliability of 
binary networks. The work by Dohmen, mentioned in the previous exam- 
ple, provides reliability bounds and the exact formula for a wider class of 
networks. However, there are other binary structures that are not networks. 
A fc-out-of-n system is well studied but it is only one of the many possi- 
ble examples of a nonnetwork system. In our setting a /c-out-of-n system 
can be represented, as any other coherent system, by a monomial ideal. For 
instance, the ideal 

M = (xiX2, 2:1X3, j;iX4, X2Xs,X2X4, X3X4) 

is the monomial ideal corresponding to the nonfailure set of a 2-out-of-4 sys- 
tem. In general, different situations can arise where the minimal nonfailure 
points do not correspond to the paths of any binary networks. For example, it 

can be seen that the minimal nonfailure set T* = {(1, 1, 0, 0, 0), (0, 1,1,0,0), (0, 0, 1, 1, 0), (1, 0,0,1,1), (0, 1 
cannot derive from any network. As in the previous example, labeling the 
minimal nonfailure points (1, . . . ,5), after deformation and ranking the re- 
sulting Scarf complex is obtained: 

A = {{1235}, {123}, {135}, {125}, {235}, {345}, 

{13}, {12}, {15}, {34}, {23}, {25}, {35}, {45}, {1}, {2}, {3}, {4}, {5}}. 

5.4. The multistate case. As explained, the full power of the method is 
exhibited for multistate systems where the literature on identities is sparser 
than for the binary case and bounds even more so. The same steps (i)- 
(iv) are used as in the binary case. For multistate systems the degree of 
nongenericity tends to be relatively lower because the multiple levels tend 
to separate out the points better. In the next section we show that with 
some care the generic setting can also be extended to cover the continuous 
case. 

In Section 5 we introduced the notion of a system without imposing any 
limit on the number of state levels that each component can assume. The 
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orthants Qa as defined above were infinite portions of the state space. How- 
ever, in the situation where each component has only a finite number of 
states the state space "D is finite and the nonfailure set JF is given by the 
intersection of the monomials in the monomial ideal generated by the set 
T* with the state space "D. The formula in (8) is still valid if we take care 
to intersect each orthant with the set T>. Full details are given in [11], where 
the finite case is described in a more formal setting. 

We consider here a multistate system where the nonfailure set is defined 
by setting a cutoff point c for an increasing function ^ defined on the state 
space V, 

f'c = {a£V\'^{a)>c}. 

An example of such a setting is given by a situation where the profit ^ 
deriving from a certain system is a function of the performance of each single 
component. In addition, we consider the case where the performance (profit) 
of the system is mostly affected by the performance of certain components 
or combination of components. For example, consider the following profit ^ 
for a given four-component system: 

(11) ^ = ai -|- 02 + 4q;3 -|- 5a4 + 2a3Q;4, 

where the performance Oi of the ith component can take values in {0, 1, 2, 3}. 
The coefficient 4 for the term 03 means that the contribution to the profit 
^ from the third component is 4 when the component is at level 1, 8 when 
the component is at level 2 and so on. The "interaction" term 20304 implies 
that profit is further boosted when components 3 and 4 are simultaneously 
at high performance levels. 

For the profit function (11) defined on the state space T> = {0,1,2,3}^, 
if the cutoff point is set at at c = 28 we obtain nine minimal nonfailure 
configurations, 

J-2*8 = {(3, 2, 3,1), (2, 3, 3,1), (2, 0,2, 2), (1,1, 2, 2), 

(0,2, 2,2), (3,0, 1,3), (2, 1,1,3), (1,2, 1,3), (0,3, 1,3)}. 

Since the profit function is increasing in all the components, any con- 
figuration P with P ^ a, with a in ^28' ^^^^^ belongs to the nonfailure set 
jr2g = {Q,e2)!^(o) >28}. 

To obtain the reliability 7^ = Prob(.F28) a-nd the corresponding bounds 
we use the improved inclusion-exclusion formula via the Scarf complex, as 
described above. Since the genericity condition does not hold for jr|g we 
need to proceed via deformation and ranking: 

J-2*8 = {(7, 4, 7,0), (4, 7, 8,1), (5, 0,4, 2), (2, 2, 5, 3), 

(0,5, 6,4), (8, 1,0, 5), (6,3, 1,6), (3,6,2, 7), (1,8,3, 8)}. 
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After labeling the minimal nonfailure points in jF|g from 1 to 9 the facets 
of the corresponding Scarf complex can be obtained: 

A28 = {{123}, {489}, {459}, {234}, {378}, {348}, {367}, 

{49}, {59}, {23}, {24}, {38}, {89}, {78}, {48}, {36}, 

{67}, {45}, {34}, {13}, {37}, {12}, {9}, {8}, {7}, {6}, {5}, {4}, {3}, {2}, {1}} . 

The geometrical representation of the Scarf complex is given in Figure 2. 
Notice that the inclusion-exclusion formula now includes only 31 terms com- 
pared to the 2^ — 1 = 511 terms that would appear in the complete classical 
formula. 

As an example of the calculation of the reliability formula (8) or the tube 
bounds (9), we obtain the probability of the orthant associated with the 
index set {489} in the Scarf complex A28- 

Using the notation introduced in [9], we indicate with Xi the random state 
of component i and with Pi{j) the probability that component i is in a state 
level "better" than j: Pi{j) = Prob(Xj > j). In addition we indicate with ptj 
the probability that component i is exactly in state j: pij = Prob(Xj = j), 
i = 1, . . . ,4, j = 0, . . . ,3. Therefore we obtain the following expression for 

my- 



n 



J 
k=j 

where rij is the maximum performance level of the j'th component. The index 
set / = {489} in the Scarf complex corresponds to the least common multiple 
of the three monomials in the (original, before perturbation) minimal failure 
sets F2S labeled 4, 8 and 9, respectively: 

TTlj ^ \.C1T1{X-\XoXoXaj X^XoX-dXa^ X^XoXa) —— X^XoXoXa- 




5 

Fig. 2. Geometric realization of the Scarf complex in the multistate example. 
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Thus, in the case of mutual independence among the components, the prob- 
abihty of the orthant associated with the index set I = {489} [and the cor- 
responding monomial mj = x^xlxlx^, a/ = (1,3,2,3)] is obtained by 

4 

Prob(Q„,) = n Prob(X, > a«) 

i=l 

3 3 

= Pi(l)^2(3)P3(2)P4(3) =P23 -^43 ' Y.P^J ' E^S,- 

5.5. The continuous distribution case. A natural generalization of the 
multistate system is to components with continuous states. With care, the 
continuous case can be mapped into the discrete case in such a way that the 
inclusion-exclusion lemma and associated tube bounds can be inferred back 
from the discrete to the continuous case. 

A simple way to define the mapping is to replace continuous values by 
their rank along each dimension. In addition one needs to "quantize" the 
continuous distributions. 

Thus let Z = (Zi, . . . , Zn) be an n-vector random variable with distribu- 
tion function 





Fz{z) = 


:Prob(Z<^). 


Define critical points 






(12) 




, . . . , 2 


and the corresponding 


orthants 






Q-A^) = 


:{Z|Z>Z«}. 


Thus the nonfailure event is given 


by 




T^ 


m 



Now consider for a particular dimension i the set of all values on that di- 
mension, 

Jl) Ira) 

^i ) • • • ! ^i 

Assume the continuous genericity (general position) condition, 

(13) for all 3 / 3 =^ zf^ + zf'^ j = 1, . . . , m. 

The condition is not too strong in this case because for continuous variables 
the critical points (12) can be easily chosen so that condition (13) is satisfied. 
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Fix the dimension i and relabel the values so that 

(1) / (2) / / (m) 

Finally, create m + 1 states for the discrete variable Xi so that 

Xi is in state j ^^ zp'^ <Zi< zp'+^^ . 
This induces a distribution on the {m + 1)" grid with probabilities 

n 
i=l 

which allows the orthant probabilities to be expressed in a straightforward 
way in terms of the distribution function F, 

Prob(Q,w) = l-F(4*\...,4)). 

The Scarf complex and associated tube bounds as obtained in the previous 
sections can then be derived. 
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